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m 

We propose an analytical approximation for the dependence of the effective charge on the bare 
charge for spherical and cylindrical macro-ions as a function of the size of the colloid and salt 
O^l . content, for the situation of a unique colloid immersed in a sea of electrolyte (where the definition of 

an effective charge is non ambiguous). Our approach is based on the Poisson-Boltzmann (PB) mean- 
er | field theory. Mathematically speaking, our estimate is asymptotically exact in the limit «a > 1, 
where a is the radius of the colloid and k the inverse screening length. In practice, a careful 
' comparison with effective charges parameters obtained by numerically solving the full non-linear 
PB theory proves that it is good down to ko ~ 1. This is precisely the limit appropriate to treat 
colloidal suspensions. A particular emphasis is put on the range of parameters suitable to describe 
^l-H both single and double strand DNA molecules under physiological conditions. 

O ' 

53 I I. INTRODUCTION 

We know from the work of Debye and Hiickcl that elementary charges, e, immersed in an electrolyte solution interact 
through a screened Coulomb pair potential: V(r) ~ e 2 exp(— nr)/r, where the screening length, k^ 1 , characterizes 
the thermodynamics of the ionic species. For highly charged macro-ions (bare charge Z^ aic e, where Zh arc 3> 1), the 
strong electrostatic coupling between the macro-ion and the micro-ions results in an additional screening so that the 
actual Coulomb pair potential still writes V(r) ~ Z^ s e 2 exp{—n'r)/r at large distances (with possibly n' ^ n [1]), but 
now Z e g is an effective (also called "apparent", or " renormalized" ) charge parameter, much smaller in absolute value 
than .Zbare- The idea is that the charged colloid retains captive a fraction of the oppositely charged micro- ions in its 
I , immediate vicinity, and therefore apparently behaves as a new entity of lower electric charge. 
\Q ' Because the effective charge is the relevant parameter to compute the electrostatics of the system at large inter- 
particle distances, the concept of charge renormalization plays a central role in the thermodynamics of highly charged 
t-H i colloidal suspensions. Several reviews have appeared recently which discuss this notion [1-4], see also [5]. In the 
colloid science field this concept has been introduced by Alexander et al. [6] in the context of the Poisson-Boltzmann 
(PB) cell model, but had been widely accepted since the fifties in the field of linear polyelectrolytes [7,8]. For an 
isolated macro-ion on the other hand, the definition of an effective charge from the far field potential created in an 
electrolyte is unambiguous [2,9-11]. 

In general the effective charge depends on the geometry of the particle, the concentration of macro-ions and the 
thermodynamics of the electrolyte. Within the non-linear Poisson-Boltzmann (PB) mean-field theory: Z e s ~ Z^axe 
"-^J ■ for low values of Z^ arc , and Z e g saturates to a constant, Z^, when Zb are — ► oo. Many studies have focused on 
finding an approximation for Z^g. In particular, in the limiting case of infinite dilution, we have proposed a matching 
procedure for the electrostatic potential such that may be estimated at the level of the linearized PB theory only 
(the non linearity of the problem occurring only through an effective boundary condition, see [9,10]). This approach 
may be generalized to account for concentrated suspensions. In contrast, much less is known about the functional 
dependence of Z e g on .Zbare- This is because finding the exact analytical dependence Z e R(Z\, are ) involves solving the 
full non-linear PB system of equations, which is out of reach in general. 

Here, we propose an analytical approximation for the dependence of Z c g on Zb ar c for spherical and cylindrical 
macro-ions as a function of the size of the colloid and the salt content. We restrict ourselves to the infinite dilution 
limit where an exact analytical representation of the electrostatic potential, ip, solution of the non-linear PB theory 
has been recently obtained [12]. Our estimate is asymptotically exact in the limiting case when > 1. In practice 
however, it is an accurate approximation of the exact solution in the whole colloidal domain: na > 1, where a is the 
typical size of the colloid. 
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II. EFFECTIVE CHARGE DEPENDENCE 



We consider the situation of an isolated macro-ion of given bare charge in a symmetrical, mono-valent electrolyte of 
bulk density no (no confinement). The solvent is considered as a medium of uniform dielectric (CGS) permittivity e. 
Within Poisson-Boltzmann theory, micro-ions/micro-ions correlations are discarded and the potential of mean force 
identified with the electrostatic potential ip. Accordingly, the reduced electrostatic potential (0 = eip/kT), assumed 
to vanish far from macro-ion, obeys the equation 



where the screening factor n is defined as n 2 — 8n£Bn and the Bjerrum length quantifies the strength of electrostatic 
coupling: £b = e 2 /(ekT) (kT is the thermal energy). Far away from the colloid, (j> obeys the linearized PB (LPB) 
equation V 2 = K 2 <fi. The far field within PB is thus the same as that found within LPB, provided the charge is 
suitably renormalized is this latter case. In particular, for infinite rods (radius a, bare line charge density Abar C e), we 
have 



where Kq is the zero order modified Bessel function (the symbol <~ stands for "asymptotically equal" ). In Eq. (2), 
the prefactor Aoq is a function of both the thermodynamics of the electrolyte and the characteristics of the macro-ion: 
^oo(Abarc, ko). With these notations, the effective line charge density, A ff, is such that 



where K\ is the first order modified Bessel function. 

Using the method of multiple scales, Shkel et al were able to propose an approximate expression for Aoo(Abarei wa) 
up to the second order in (na)^ 1 <C 1 [12]. This result can be translated into an approximate form for the functional 
dependence of A e ff(Abare) through Eq. (3). After some algebra, we find, again for cylinders and up to order 0(l/(na)) 



V 2 (f> = k 2 sinh0, 
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where 
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and the function T, also useful in spherical geometry (see below), is defined as 
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In the limit of diverging bare charge Ab; 
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1 so that the saturation value for the line charge density reads 




(7) 



In the opposite limit where Abarc — * 0, Eq. (4) yields A e ff = Abarc as expected. 
Similarly, for spheres (radius a, bare charge i?ba re e) we find that 




(8) 



where now 
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and the function T{x) is still defined by Eq. (6). The corresponding saturation value for the effective charge is 



1 



while for low bare charges, we have 



Zg=4Ka + 6 + 0[ — ), (10) 



lim -^L = 1. (11) 

■^barc^O ^barc 



Expressions (4)-(7) [resp. Eq. (8)-(10)] provide the first analytical estimate of the functional dependence of the 
effective charge on the bare charge for cylindrical (resp. spherical) macro-ions in an electrolyte solution. They are 
the exact expansions of the correct result in the limit of large na. In practice however, they are accurate as soon as 
na > 1, as we now show. 

We compare on Figs. 1 and 2 the results for the effective charge against the bare charge obtained using: a) our 
estimate and b) the exact solution found by solving the full numerical PB problem. The results are displayed for two 
different values of the ionic strength of the solution. The agreement is seen to be excellent for na > 5, and reasonable 
down toKoxl (see also below) . 



III. CYLINDRICAL GEOMETRY AND DNA 



Because DNA (a rod-like polyelectrolyte in first approximation) is of paramount importance in biology, the cylin- 
drical case deserves a special attention. For typical parameters of double strand DNA (a = 10 A, Abarc^B — 4.2), the 
domain na > 1, where we expect our estimate to be good enough, reads no > 0.1 M for monovalent symmetrical 
electrolyte. This is experimentally relevant. Indeed, physiological conditions are found for n = 0.15 M (k" 1 = 8 A). 
For simple strand DNA (a = 7 A, Abarc^s — 2.1) the condition na > 1 writes no > 0.2 M. Fig. 3 displays the effective 
charge, Aeff^B, against na for Abaro^s = 2.1, and Abaro^s = 4.2. We observe that our estimate remains decent down 
to na — 0.5, where the error is of the order of 7% (see also Tables 1 and 2). 

Note that, for DNA in typical experimental conditions, the parameters are such that we find ourselves in the 
crossover regime where both the linear approximation (A e ff = Abare) and the asymptotic approximation (A er i = A||j) 
fail to be accurate, as shown in Tables 1 and 2. In other words, both Abare and A|fj provide a rather poor approximation 
for Aefr- This justifies a posteriori our effort to find a good approximation for the full functional dependence of 

Acff( Abare)- 



IV. CONCLUDING REMARKS 



In this contribution, we have proposed an analytical approximation for the dependence of the effective charge on 
the bare charge for spherical and cylindrical macro- ions as a function of the size of the colloid and salt content. 
Mathematically speaking, our estimate is asymptotically exact in the limit na >• 1. In practice, a careful comparison 
with effective charges parameters obtained by numerically solving the full non-linear PB theory proves that it is good 
down to na > 1 (where a is the radius macro- ion). This is precisely the relevant range of parameters for colloidal 
suspensions. 

An important example considered in some detail is the cylindrical geometry This is because the infinitely long 
charged rod provides a simple model for DNA. While the simple asymptotical approximations fail for DNA in typical 
physiological conditions, we provide an estimate for the effective charge line density, which compares well with the 
numerical results. 

We have performed the analysis at the level of the mean-field Poisson-Boltzmann theory. In spite of its limitations, 
this picture is excellent for all existing macro-ions in water at room temperature, when only monovalent micro-ions are 
present in the electrolyte (see e.g. ref [4] for an estimation of a relevant coupling constant quantifying the importance 
of the micro-ionic correlations). It may also hold for multivalent counter- ions, provided the surface charge density of 
the macro-ion is not too large. In general, for a given total charge of the macro-ion, micro-ionic correlation become 
irrelevant in the asymptotic limit where n» £b [13,14]. 
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FIG. 1. Rescaled effective charge, Z^ts/a, versus rescaled bare charge, ^bare^s/o, for isolated spherical macro-ions in a 
symmetrical, monovalent electrolyte solution (for na — 5, where k" 1 is the Debye-Hiickel screening length of the electrolyte, 
and a is the sphere radius). The solid line is the analytical estimate of Z e f{ [see Eq. (8)], whereas the open circles are the exact 
values of Z c g found by numerically solving the full non-linear PB problem. The dashed line has a of slope unity to emphasize 
the initial linear behaviour. The inset shows ^off(^barc) f° r Ka = 1- 
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FIG. 2. Effective line charge density, X c g, versus bare line charge density, A^aro f° r an isolated, infinitely long cylindrical 
macro-ion (/to = 5). The solid line is the analytical estimate of A off [see Eq. (4)], whereas the open circles are the exact values 
of A e ff found by numerically solving the non-linear PB theory (the dashed line has a of slope unity). The inset shows A e ff(Abare) 
for kcl — 2. 
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FIG. 3. Effective line charge density, A e g-£s, versus na for isolated, infinitely long cylindrical macro-ions of bare charge 
Abare^s = 2.1 (lower sets) and A barc ^B = 4.2 (upper sets). The solid line is the analytical estimate of \ c g [see Eq. (4)], whereas 
the open circles are the exact values of A e g found by numerically solving the full non-linear PB theory. In the formal limit 
ko — ► oo, these two curves would saturate to 2.1 and 4.2. 
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Table 1: Comparison of the effective charge and the saturation value for single strand DNA at physiological conditions 
as found a) by our analytical estimate, b) numerically solving PB theory 



Abarc^B 


KCI 


X cS £ B [Eq.{4)\ 


Acff^B (num) 


Xgis [Eq.(7)\ 


\ s e fftB (num) 


2.1 


0.25 


1.35 


1.20 


2.00 


1.61 


2.1 


0.5 


1.50 


1.40 


2.50 


2.20 


2.1 


1 


1.69 


1.65 


3.50 


3.27 


2.1 


1.5 


1.81 


1.78 


4.50 


4.32 


2.1 


2 


1.88 


1.87 


5.50 


5.35 



Table 2: Same as Table 1 for double strand DNA. 



Abarc^B 


KCI 


X cS £ B [Eq.{4)\ 


Aoff^B (num) 


Xge B [Eq.(7)\ 


X s e $t B (num) 


4.2 


0.25 


1.65 


1.41 


2.00 


1.61 


4.2 


0.5 


1.94 


1.78 


2.50 


2.20 


4.2 


1 


2.40 


2.31 


3.50 


3.27 


4.2 


1.5 


2.76 


2.71 


4.50 


4.32 


4.2 


2 


3.02 


2.99 


5.50 


5.35 
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